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LEVEL  SHIFTS  AND  INELASTIC  ELECTRON  SCATTERING  IN  DENSE  PLASMAS 


I.  INTRODUCTION 

The  plasma  polarization  shift  of  spectral  lines  remains  a  basic  and 
controversial  effect.  Grlem^  tried  to  explain  the  blue  shifts  for  the 
resonance  lines  of  ionized  helium  by  considering  the  perturbation  of  the 
upper  levels  by  free  plasma  electrons.  An  alternate  approach  Involves  the 

II 

solution  of  the  Schrodinger  equation  for  the  bound  electrons  in  the  Debye- 

It 

Huckel  potential.  This  leads  to  red  shifts  of  the  lines  because  it  causes 

large  shifts  of  the  lower  levels.  Experimentally,  evidence  has  been  reported 
2  3 

for  blue  shifts  of  Lyman  lines  and  red  shifts  for  Balmer  lines  for  He  II, 
but  shifts  measured  for  other  similar  ionized  emitters  are  much  smaller  or 
zero. 

4  " 

Theimer  and  Kepple  have  shown  that  the  theory  using  the  Debye-Huckel 

potential  is  inadequate  because  it  ignores  the  reaction  of  plasma  electrons 
to  the  presence  of  the  bound  electron.  Their  self-consistent  calculation 
for  a  hydrogen  plasma  indicates  much  smaller  level  shifts  than  those  pre¬ 
dicted  by  previous  theories.  More  recently,  an  analogous  treatment  was 

5  6 

adopted  by  Rozsnyai  and  Skupsky  with  similar  results.  In  these  theories 
the  free  electron  gas  was  treated  classically  and  the  energy  levels  were 

II 

obtained  from  the  solution  of  the  Schrodinger  equation. 

In  the  present  investigation  we  obtained  level  shifts,  transition  prob- 

*4-9  4*17 

ability  coefficients  and  collision  cross  sections  for  He  and  Ar  in  a 

24  24  -3  +9 

plasma  with  electron  density  ranging  from  10  to  6  x  10  cm  (for  Ne  ) 

25  25  -3  +17 

and  from  2  x  10  to  8  x  10  cm  (for  Ar  )  using  a  quantum-mechanical 
description  of  the  plasma  electrons.  The  theory  in  a  simplified  form  was 
applied  to  the  low-density  region  and  level  shifts  were  obtained  for  H  and  He+. 


II.  DESCRIPTION  OF  THE  METHOD 

We  consider  an  isolated  impurity  ion  with  nuclear  charge  Z  immersed 
in  a  fully  ionized  hydrogen  plasma.  For  simplicity  all  electric  charges  are 
assumed  to  have  a  spherically  synmetrlc  charge  distribution.  In  addition, 
the  positive  charge  due  to  protons  is  continuously  distributed.  This  last 
assumption  is  probably  the  most  drastic  simplification  that  we  have  adopted. 

Manuscript  submitted  October  6, 1981. 
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Two  types  of  ions  will  be  studied:  completely  stripped  ions  end  one- 
electron  hydrogenic  ions.  In  the  letter  case,  it  is  es sumed  that  in  a  dense 
hot  plasma  the  average  local  charge  density  due  to  the  bound  electron  can 
be  written  as  (atomic  units  are  used  throughout  the  paper,  unless  indicated) 
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where  b  are  the  occupation  numbers,  Pnjl(r)  are  normalized  radial  wave 
functions,  and  b  ^  *1.  In  order  for  this  assumption  to  be  valid  it  is 

n  ii 

necessary  that  the  average  time  between  collisions  leading  to  ionization  and 
recombination  be  substantially  larger  than  the  time  between  collisions  that 
change  .he  state  of  excitation.  Generally,  the  occupation  numbers  represent 
an  average  state  of  excitation  of  the  ion  and  are  determined  by  all  possible 
colllsional  and  radiative  processes  involving  the  bound  electron.  Due  to 
the  complexity  of  this  situation,  simplifying  assumptions  about  b  will  be 
adopted. 

If  p  and  p  represent  local  charge  densities  due  to  the  free  electrons 
e  p 

and  protons,  the  electrostatic  potential  V(r)  at  the  distance  r  from  the 
nucleus  with  charge  Z  is  then  given  by 


V(r) 
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At  large  distances  from  the  ion,  •  0  and  the  plasma  is  neutral,  so  that 


P 


e 


(3) 


The  positive  charge  distribution  is  described  by  classical  Boltzmann 
statistics,  viz. 

Pp  "  om  exp  (-V/ltgT) ,  (4) 

where  V(r)  is  given  by  (2) ,  T  is  the  temperature  and  the  Boltzmann 
constant. 
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The  free  electron!  ere  represented  by  e  Fertui-Direc  energy  dlstrlbulton 

outside  e  spherical  boundary  corresponding  to  the  distance  r  from  the 

o 

nucleus;  hovever,  Inside  this  boundary  they  can  be  treated  quantum-mechani- 
cally  and  described  by  wave  functions  which  are  solutions  of  the  time  in- 

II 

dependent  Schrodinger  equation.  In  practical  calculations  the  distance  r  is 

o 

taken  to  be  large  enough  so  that  plasma  at  the  boundary  may  be  considered 

neutral  and  condition  (3)  valid  within  numerical  accuracy.  The  wave  functions 

1  2 

for  the  free  electrons  with  energy  jk  at  r  >  tq  can  be  expanded  in  partial 
waves  and  the  expression  for  the  spherically  averaged  charge  density  due  to 
the  free  electrons  at  any  r  is  then  given  by 
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-  J  W(k)  (kr2)”1 
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(21+1)  F, 


kJ> 


(r)  dk. 
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where 


W(k)  -  7T  ~2k2 


1  +  exp 


(6) 


y  is  the  chemical  potential  of  the  free  electron  gas  determined  by  the 
condition 

00 

J  W(k)  dk  - 

o 

and  the  functions  F^  are  solutions  of  the  equation 


(7) 
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+  2V(r)  -  2V  (r)  -  2V _ (r) 


dr 


ex 


corr 


(8) 


+  2V  (»)  +  2V  (®)  +  k*  F.  „(r)  -  0 
ex  corr  J  k£ 


with  the  asymptotic  form 

_  1 

FkJl(r)  a,  k”  2  sin  (kr  +  <Sk£) .  (9) 
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The  electrostatic  potential  V  la  given  by  (2)  and  the  exchange  energy  V  and 

OX 

the  correlation  energy  V  of  a  free  electron  in  a  plasma  vara  taken  from 

7  8  corrg 

Dharma-vardana  '  and  Taylor  and  are  given  in  terms  of  temperature  T  and 

local  electron  density  N  =  -  n  : 

e  e 

V  -  -  0.4073  r  -1  tanh  (t-1) ,  (10) 

OX  8 

%i 

V  -  -  0.6109  r  2  (-0.0081  +  1.127 t 2  +  2.756x4) 

corr  s 

1  (11) 
x(1.0  +  1.291t2  +  3.593x4)”^  tanh  (x  '*), 


x  -  2  ( — ^ — )  3 
9tt 
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In  practice,  the  Integral  in  (5)  is  replaced  by  a  summation  and  an  8-point 

Gaussian  quadrature  was  used  to  evaluate  p  for  each  value  of  r.  At  large 

distances,  the  convergence  of  the  sum  over  partial  waves  in  (5)  is  rather 

slow  and  we  found  it  necessary  to  include  up  to  50  or  70  partial  waves, 

depending  on  the  value  of  k,  to  achieve  sufficient  accuracy  for  r  <  r  . 

o 

The  effect  of  Fermi  degeneracy  is  taken  Into  account  by  the  distribution 
function  W(k)  in  (5)  for  densities  equal  to  p^.  For  the  cases  studied  in  this 
paper  the  Fermi  degeneracy  causes  only  minor  deviations  from  the  Maxwell- 
Bolt  zmann  dlstribuiton .  Inside  the  boundary  rQ,  where  the  plasma  electrons 
in  our  model  are  described  by  wave  functions  and  where  the  local  electron  den¬ 
sities  may  reach  values  substantially  larger  than  p  ,  the  effect  of  degeneracy 
should  be  properly  taken  into  account  by  using  an  antisymmetric  function  for 
tne  whole  system.  Due  to  complications  arising  from  this  procedure,  we  have 
ignored  the  possible  additional  effects  of  Fermi  degeneracy  inside  rQ.  Also 
all  effects  due  to  relativistic  corrections  have  been  ignored  since  they  are 
probably  negligible  for  light  nuclei. 

Sadial  functions  of  the  bound  electron  sacisfy  the  equation 


2V,  (r)  +  2E  . 
b  n,. 


F  ,(r>  •  0. 


I 


(12) 


The  potential  acting  on  the  bound  electron  is  given  by  (2)  with  P ^  *  0. 

The  energy  shift  of  a  hydrogenic  level  nt  due  to  the  plasma  polarisation  is 

1  En*  •  En*  +  22(2“2>‘1  '  <“> 

Equations  (1) ,  (2) ,  (4) ,  (5) ,  (8)  and  (12)  should  be  solved  self-con- 
sistently  with  the  boundary  condition 

rV(r>  -+  0,  PnA(r)  0  (14) 

for  r  -*■  «.  In  practice  this  condition  should  be  satisfied  at  r  ■  r  . 

The  equations  were  solved  iteratively  starting  with  hydrogenic  P  .  (for 

i»  n* 

hydrogenic  ions)  and  with  V(r)  Initially  replaced  by  a  Debye-Huckel  potential. 
First  we  found  a  self-consistent  potential  (2)  while  keeping  P  „  unchanged. 
This  procedure  requires  up  to  four  iterations.  Then  we  used  the  new  form  of 
P  ^  obtained  from  (12)  for  the  following  sequence  of  iterations.  It  was 
found  that  in  all  cases  presently  studied  the  second  sequence  of  iterations 
did  not  lead  to  any  substantial  improvement  of  results  and  that  the  first 
sequence  already  yielded  self-consistent  values  of  P^  and  E^. 


HI.  RESULTS  FOR  Nc*10,  Ne*9  and  Ar*17 

Figure  1  shows  two  typical  examples  of  the  quantity  rV  for  Ne+i° ,  where 
V  is  the  self-consistent  potential  (2)  with  p.  »  0.  These  results  are  com- 

if  0 

pared  with  the  Debye-Huckel  theory  of  plasma  screening,  which  agrees  well 
with  the  self-consistent  potential  at  distances  larger  than  2D,  where  D  is 

it 

the  Debye  screening  length,  but  the  Debye-Huckel  potential  is  systematically 
lower  in  the  inner  region.  The  difference  appears  to  increase  with 
decreasing  values  of  D. 

The  calculations  for  hydrogen-like  Ne  and  Ar  were  performed  with 
the  simplifying  assumption  that  the  levels  Is,  2s  and  2p  are  populated 
according  to  Boltzmann  statistics  and  that  all  occupation  numbers  b  ^  with 

n  >  2  are  zero.  The  local  densities  n  (r)  and  n  (r)  of  free  electrons  and 

24-3  e  P 

protons  for  N  -  3  x  10  cm  and  T  •  200  eV  are  displayed  on  Figs.  2  and  3 
+9  * 

for  Ne  .  The  oscillatory  behavior  of  ng  is  caused  by  interference  effects 
and  is  mainly  due  to  the  slowly  moving  electrons. 
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T  -  200  eV 


r<»o> 

+0 

Figure  2.  Density  distribution  of  free  electrons  and  protons  for  Ne  . 

Densities  are  normalized  to  unity  at  infinity. 

N  -  3  x  102\  T  -  200  eV. 
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Table  1  presents  the  atomic  level  shifts  (13)  for  two  temperatures  and 
three  electron  densities  Resulting  line  shifts  of  the  Lyman  a  radiation 

are  shown  In  Figs.  4  and  S.  The  line  shifts  are  assumed  to  be  equal  to  the 
difference  of  level  shifts.  The  present  values  of  line  shifts  are  slightly 
larger  than  the  results  of  Skupsky^  who  treated  the  free  electrons  classi¬ 
cally  with  quantum  mechanical  corrections  to  the  potential. 

Table  2  shows  how  the  Einstein  coefficients  of  transition  probability 
A  for  the  Lyman  a  line  are  affected  by  plasma  conditions.  ^ 

coefficient  corresponding  to  the  pure  Coulomb  field,  i.e.  to  *  0.  The 
change  of  transition  probabilities  is  a  consequence  of  the  line  shift  combined 
with  the  change  of  radial  functions  P  .  In  order  to  ascertain  the 
sensitivity  of  line  shifts  to  the  choice  of  b  . ,  we  carried  out  additional 
calculations  for  Ne  with  the  assumption  that  only  the  Is  ground  level  is 
occupied.  The  resulting  line  shifts  were  found  to  be  only  about  1%  smaller 
than  those  shown  in  Fig.  4. 

In  the  present  calculations  we  have  ignored  all  effects  related  to  the 
time  dependent  perturbations  by  plasma  electrons  and  protons.  In  particular, 
a  possible  redistribution  of  electron  velocities  inside  the  boundary  rQ  due 
to  electron-electron  and  electron-proton  collisions  was  not  taken  into 
account  and  effects  leading  to  broadening  of  spectral  lines  were  also  omitted. 
Thus  the  results  indicate  the  magnitude  of  effects  caused  by  plasma  polari¬ 
zation  alone.  In  situations  where  there  are  only  fev.T  electrons  inside  the 
Debye  sphere  and  collisions  do  not  happen  very  often,  the  present  treatment 
should  be  perferable  to  the  procedures  based  on  statistical  considerations 
for  plasma  electrons. 

The  calculation  of  collision  excitation  cross  sections  and  rates  in  a 
dense  plasma  should  proceed,  in  principle,  along  the  lines  of  many-body 
quantum  theory.  However,  the  two-body  approximation  may  still  be  useful 
not  only  for  approximate  evaluation  of  excitation  rates,  but  especially 
because  its  simplicity  allows  us  to  separate  various  plasma  effects. 

The  change  in  the  value  of  excitation  cross  sections  and  rates  in  dense 
plasmas  is  the  result  of  several  causes:  (a)  changed  excitation  energy, 

(b)  changed  wave  functions  of  atomic  electrons,  (c)  distortion  of  the 
incident  and  scattered  waves  by  the  charge  cloud  around  the  target  ion, 

(d)  modification  of  the  Coulomb  interaction  potential  jr^  -  s  r^  1  for 

the  colliding  and  bound  electron  by  plasma  effects,  and  (3)  other  effects  of 
many  -b  ody  int  e  r  act  ion  s . 
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Table  1 


LEVEL  SHIFT  (eV) 


T(eV) 

N  (cm"3) 

e 

Is 

2s 

2p 

Ne  X 

200 

1024 

122.8 

115.6 

117.5 

3xl024 

194.0 

173.3 

178.5 

6xl024 

225.5 

189.8 

198.1 

500 

1024 

90.0 

85.5 

86.7 

3xl024 

147.3 

134.2 

137.5 

6xl024 

199.6 

174.3 

180.3 

- 

Ar  XVIII 

1000 

2xl025 

475.8 

439.9 

449.1 

5xl025 

702.0 

616.4 

637.5 

8xl025 

865.0 

736.0 

766.4 

2000 

2xl025 

375.8 

351.0 

357.3 

j 

5xl025 

563.8 

504.5 

518.8 

8::1025 

695.0 

602.6 

624.3 

10 
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Figure  5.  Density  dependence  of  the  Ar  XVIII  Lymn  -  a  line  shift 


Table  2 


A(2p  -  Is) /A.  t1  (2p  -  la) 


N  (cm  3) 

T(eV) 

e 

1024 

3  x  1024 

6  x  1024 

Ne  X 

200 

0.967 

0.904 

0.840 

500 

0.981 

0.942 

0.887 

N  (cm-3) 
e 

25 

2xlOz;> 

25 

5xK/ 

8xl025 

Ar  XVIII  1000  0,950  0.879  0.816 


2000 


0.967 


0.919 


0.873 


We  have  calculated  collision  strengths  for  the  Is  -  2s,  Is  -  2p,  and 

+17 

28  -  2p  transitions  of  Ne  and  Ar  using  the  formalism  of  the  distorted- 
wave  method  without  exchange.  Excitation  energies  and  atomic  wave  functions 
were  taken  from  our  self-consistent  calculations  and  the  equation  for  the 
radial  functions  of  the  colliding  electron  were  solved  using  the  self-con¬ 
sistent  potential  (2).  However,  the  modificaiton  of  the  potential  of  mutual 
Interaction  r^2~^  cannot  be  properly  described  within  the  formalism  of  the 
distorted-wave  approximation,  unless  the  wave  functions  of  plasma  electrons 
take  into  account  correlation  effects  between  the  bound  and  the  free  elec¬ 
trons.  The  variational  expression  for  the  eletjent  of  the  reactance  matrix 
then  contains  an  additional  term  that  partially  cancels  the  term  of  direct 
Coulomb  interaction  r^2  ^  as  a  result  of  decreased  local  density  of  plasma 
electrons  in  the  vicinity  of  the  bound  electron.  One  can  approximate  this 
effect  by  using  the  free  electron  wave  functions  without  correlation  and 
replacing  rl2_1  ***  the  matrix  ©leme^t  by  r12_1  f^rl2^’  where  ^ri2^  is  a 
properly  defined  screening  function. 


r12~J’  can  be  exPanded  in  the  ^orm 


r,  0-1  f(r,  0)  -  r”1  X) 


r12  *'r12'  "  r>  P,  (cos  w) 

>  A  ^ 


£  A,  (r  )  (r  /r  )'A 

\]i  >  <  > 

p  ■  0 


a  T  *-y 


where  w  is  the  angle  between  r^  and  r2  and  the  coefficients  A^u  depend  on  the 
form  of  f(r^2)* 

To  demonstrate  the  effect  of  screening  of  the  mutual  interaction,  we 
n 

have  adopted  a  Debye-Huckel  form  of  the  screening  function.  The  expansion 
(15)  have  the  form 

r12_1  exp(-r12/D)  -  £  r”1  +  ^  (r>/D)2  (r  2 /r^  )  +  . . .  j  exp(-r>/D)PQ(cos  w) 

r  -j  (16) 

+  [  (1  +  r>/D)  (r</r2)  +  Jo  (1  +  r>/D)  <r>/D>2  <*</*>>  +  “-J 

x  exp  (-r>/D)  (cos  u)  +  ... 

The  results  for  selected  values  of  T  and  N  are  shown  on  Figs.  6-11.  Also 
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Figure  6.  Collision  strength  ft  (la2S,  2s2S)  for  Ne*^.  T  »'  200  eV. 

1  .  H  «  0;  2,3,4  -  »e  ■  3xl024  cm-3;  2  -  D  -  ®;  3  -  D  »  A; 
4  -  D  ■  t  A  (A-  Debye  length). 
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Figure  7.  Collision  strength  ft  (ls2S,  2p2P)  for  Ne*®.  T  -  200  eV. 

1  -  Ne  -  0;  2,3,4  -  -  3xl024  cm'3;  2  -  D  -  3  -  D  -  A; 

4  -  D  ■  j  A  (A  ■  Debye  length). 
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Figure  8.  Collision  strength  ft  (2n^S,  2p“P)  for  Ne4^ .  T  -  200  eV; 

N  ■  3  x  10  H  cm  .  l  -  D  -  2  -  D  ■  Debye  length. 
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E  (keV) 

Fif/ure  9.  Collision  strength  fi  (ls2S,  2s2S)  for  Ar+32.  T  ■  7.000  eV. 
1  -  Ne  «  0;  2,^4  -  -  8xl025  cm"3;  2  -  D  -  ®; 

3  -  D  ■  A,  A  -  D  -  —A  (A  -  Debye  length). 
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E  (keY) 

2  2  +17 

Figure  10.  Collision  strength  ft  (Is  S,  2p  P)  for  Ar  .  T  -  2000  eV. 

1  -  N  -  0;  2,3,4  -  N  -  8  x  1025  cm"3;  2  -  D  «  ®; 

ft  4  6 

3-D  -  A;  4  -  D  »  j'  A  (A  ■  Debye  length) . 


19 


—  -  ^  -  -  - — I  - 


ro 


-S,  2)0 


E  (keV) 

2  2  +17 

Figure  11.  Collision  strength  ft  (2s  S,  2p  P)  for  Ar  .  T  ■  2000  eV; 

N  ■  8x10^  cm”"*.  1  -  D  ■  »;  2  -  D  ■  Debye  length. 
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shown  am  the  collision  strengths  without  any  plasma  affacts  (N#  ■  0)  for 
tha  transitions  Is  -  2s,  Is  -  2p.  The  change  of  atomic  wave  functions  by 
plasma  polarization  has  a  very  small  affect  on  the  cross  sections  for 
conditions  considered  in  this  investigation.  The  change  of  excitation 
energies  produces  a  negligible  effect  on  the  Is  -  2s  and  Is  -  2p  transitions 
but  it  removes  the  degeneracy  of  the  2s  and  2p  levels  and  reduces  the  2s  -  2p 
cross  sections  to  finite  values.  For  the  Is  -  2s  and  Is  -  2p  transitions, 
calculations  were  made  for  3  values  of  the  screening  parameter  D  in  (16) . 

D  ■  «  corresponds  to  the  situation  where  the  modification  of  mutual  inter¬ 
action  is  ignored  and  the  lowering  of  the  collision  strength  is  caused  only 
by  the  defocuslng  of  colliding  electrons  (curve  2  on  Figs.  6,  7,  9,  1C). 
Curves  3  and  4  were  obtained  with  D  ■  A  and  0  -  1/2  A,  where  A  is  the  Debye 
length  for  the  corresponding  temperature  and  density. 

The  effect  of  screening  on  partial  collision  strengths  increases  with 
increasing  angular  momentum  (Fig.  12)  and  therefore  the  dipole  transition 
Is  -  2p  is  more  affected  than  the  monopole  transition  Is  -  2s.  The  term 
corresponding  to  u  ■  1  in  the  expression  (15)  and  (16)  has  only  a  very  small 
effect  on  the  collision  strengths . 

Collisions  strengths  for  the  2s  -  2p  transition  were  calculated  for 
D  «  «>  and  D  ■  A.  The  sensitivity  of  fl  to  the  value  of  D  is  much  greater  than 
for  the  Is  -  2s  and  Is  -  2p  transitions,  because  the  excitation  energy  is 
much  smaller,  many  more  partial  waves  contribute  to  ft  at  a  given  energy  and 
higher  angular  momenta  are  more  affected  by  screening.  Collision  strengths 
exhibit  a  resonance-like  behavior  in  the  low-energy  region  above  the 
excitation  threshold.  The  enhancement  of  ft  is  caused  almost  exclusively  by 
the  contribution  from  the  l  »  2,3  incident  partial  waves  in  Ne  and  by 
2,  ■  3  in  Ar+^ , 


IV.  RESULTS  FOR  HYDROGEN  AND  IONIZED  HELIUM 

At  the  present  time ,  there  are  no  experimental  data  that  could  be  com¬ 
pared  with  the  theoretical  results  for  Ne  and  Ar+^ .  Therefore  we  have 
extended  the  calculations  to  neutral  hydrogen  end  singly  ionized  helium  which 
have  been  studied  experimentally.  Unfortunately,  the  measurement  of  plasma 
polarization  shifts  is  a  difficult  problem  and  results  of  various  experiments 
do  not  always  agree  with  each  olher.  Moreover,  the  measured  shifts  represent 
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a  combined  effect  of  plasma  polarization  (i.e.  space  charges  around  the  atom), 
ion  quadrupole  shifts,  quadratic  Stark  effect  and  shift  arising  from  electron 
collisions.  These  additional  effects  have  been  ignored  in  our  present  investi¬ 
gation  and  therefore  our  results  for  H  and  He+  have  only  a  qualitative 
character. 

The  plasma  conditions  in  which  neutral  hydrogen  and  ionized  helium  radiate 
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correspond  to  such  low  electron  densities  (£  10  cm  )  that  several  sub¬ 
stantial  slmplications  in  our  theoretical  treatment  are  possible.  In  the 
limit  of  low  charge  density  (due  to  plasma  electrons  and  positive  ions  in 
the  vicinity  of  the  radiating  atom)  we  may  omit  the  effect  of  these  space 
charges  on  the  motion  of  free  electrons  and  we  solve  eq.  (8)  with 

V  «  V  *  0  and  with  V  from  (2)  with  p  *  p  ■  0.  p  .  is  again  given  by 

ex  corr  e  p  b  ° 

(1) ,  but  due  to  the  low  temperature,  the  population  of  excited  levels  is  so 
small  that  it  may  be  omitted  altogether. 

We  will  consider  pure  hydrogen  and  pure  helium  plasma  ana  for  the  sake 
of  simplicity  we  will  assume  that  the  fraction  of  doubly  ionired  helium  ions 
is  very  small.  Then  all  positive  ions  are  singly  ionized  and  we  may  replace 
Pp  by  p.^  in  eq.  (2)  -  (4).  For  the  positive  charge  distribution  we  again 
adopt  eq.  (4)  and  in  order  to  demonstrate  the  effect  of  different  ion  distri¬ 
butions  on  level  shifts  we  use  three  different  forms  of  V:  (a)  V  *  0  which 
corresponds  to  a  uniform  ion  distribution,  (b)  V  calculated  from  (2)  with 
pe  *  Pp  *  0  311(1  with  pfc  from  (1)  assuming  that  only  the  Is  level  is  occupied, 
(c)  V  from  (2)  with  pe  "  Pp  m  0  311(1  from  (1)  assuming  that  only  the  level 
n£,  for  which  the  shift  is  being  calculated,  is  occupied.  This  last  choice 
may  find  a  justification  if  we  compart  de  Broglie  wavelength  Xg  for  thermal 
electrons  and  ions :  while  Xg  for  electrons  is  comparable  or  larger  than  the 
dimensions  of  the  atom  and  therefore  plasma  electrons  do  not  adjust  to  in¬ 
dividual  atomic  states  nil,  Xg  for  ions  is  much  smaller  and  the  ion  distribution 
may  depend  on  nJl.  In  using  eq.  (4)  for  the  ion  distribution,  one  should 
employ,  of  course,  the  true  interatomic  potential  V.  This  potential  may  be 
substantially  different,  especially  for  hydrogen,  from  the  potential  we  are 
using  the  present  work.  This  is  another  reason  why  out  results  should  be 
compared  with  experiment  only  with  caution. 

The  level  shifts  in  the  low-density  limit  are  calculated  from  the 
perturbation  expression 


(17) 


J  En .1  *  -  /  ?nt  V'  dr- 
o 

where  V*  ■  V-Zr  ^  and  V  is  calculated  from  (2)  with  p.  ■  0,  p  ■  p.  and  with 

D  pi 

and  pg  obtained  by  the  procedure  just  described.  In  the  low-density  limit 
we  do  not  perform  self-consistent  calculations  and  therefore  we  do  not  know 
the  correct  behavior  of  V  at  large  distances.  Consequently  ve  cannot  carry 
out  the  integration  in  (17)  up  to  the  indicated  upper  limit  and  cannot 
calculate  absolute  level  shifts  as  In  the  high  density  case.  However,  we  can 
obtain  relative  shifts  -  A  E^,^,  if  we  integrate,  in  (17),  up  to  the 

distance  beyond  which  all  radial  functions  P  ^  involved  in  the  calculation 
vanish . 

From  (17)  it  follows  that  level  shifts  are  proportional  to  electron 

density  Ng.  As  a  consequence,  we  do  not  have  to  perform  a  thermal  average 

of  electron  density  distribution  before  calculating  the  potential  from  (2) 

and  shifts  from  (17) ,  but  instead  we  can  first  calculate  shifts  for  the 

charge  distribution  corresponding  to  a  given  velocity  of  electrons  and  perform 

the  thermal  average  of  these  shifts.  They  are  a  smooth  function  of  electron 

velocity  and  the  average  can  be  obtained  with  great  accuracy. 

We  have  obtained  relative  level  shifts  for  all  nil  levels  of  H  and  He 

with  n  -  1  to  4.  Weighted  averages  over  the  £  -  values  for  each  principal 
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quantum  number  n,  normalized  to  the  density  »  10  cm  are  displayed  on 
Figs.  13  and  14.  Relative  level  shifts  for  H  are  more  sensitive  to  the  choice 
of  the  ion  distribution  than  for  He+. 

For  He+,  our  results  predict  a  very  small  red  shift  of  Lyman  series  lines. 

2 

The  experimental  evidence,  however,  points  to  a  much  larger  blue  shift  . 

9  10 

Theoretical  blue  shifts  were  derived  by  Greig  et  al.  ,  Burgess  and  Peacock 
and  Volonte^,  but  the  apparent  agreement  of  all  these  theoretical  results 
with  experiment  is  basically  due  to  the  omission  of  space  charges  outside 
the  orbit  of  the  bound  electron,  and  this  assumption  does  not  appear  to  be 
properly  justified. 

+ 

For  the  3-2  transition  of  He  at  1640  A  our  theory  predicts  a 

negligible  red  shift  in  agreement  with  experiment,  but  the  theoretical  shift 

o 

for  the  4-2  transition  at  1215  A  is  again  too  small,  except  possibly  for 
the  temperature  around  4  eV. 
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Figure  13.  Temperature  dependence  of  relative  level  shifts  for  hydrogen. 

Dotted  lines:  uniform  ion  distribution  (ion  density  obtained 
form  eq.  (4)  with  V»0) .  Broken  lines:  ion  density  given  by 
eq.  (4)  with  V  derived  from  screening  by  the  ni.  electron. 
Solid  line:  ion  density  given  by  eq.  (4)  with  V  derived  from 
screening  by  the  Is  electron. 
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Figure  14.  Temperature  dependence  of  relative  level  shifts  for  helium 

Dotted  lines:  uniform  ion  distribution  (ion  density  obtained 
form  eq.  (4)  with  V-0) .  Broken  lines:  ion  density  given  by 
eq.  (4)  with  V  derived  from  screening  by  the  nS.  electron. 
Solid  line:  ion  density  given  by  eq.  (4)  with  V  derived  from 
screening  by  the  Is  electron. 


^2 

Pittman  et  al.  measured  a  red  shift,  of  0.17  %  for  the  transition 

4  -  3  at  4686  &  for  T  ■  4  eV  .  From  Fig.  14  we  find  shifts  of  0.07, 

o 

1.12  and  0.14  A  for  three  different  ion  distributions.  This  may  be  con¬ 
sidered  a  satisfactory  agreement  with  respect  to  the  uncertainty  of  the  ion 
distribution. 

The  same  authors  also  measured  red  shifts  of  the  Balmer  series  lines  cf 
neutral  hydrogen.  Our  theoretical  results  predict  either  much  an-al1  values 
or  even  blue  shifts,  but  an  agreement  can  hardly  be  expected  due  to  cne  very 
approximate  form  of  the  ion  distribution. 

For  reasons  mentioned  earlier,  the  comparison  with  experiment  is 
generally  very  inconclusive.  The  theoretical  treatment  should  be  improved 
first  of  all  by  using  correct  ion  distribution  and  by  properly  taking  into 
account  other  effects  leading  to  the  shifts  of  atomic  levels. 


ACKNOWLEDGMENT 


This  work  was  supported,  in  part,  by  the  Office  of  Naval 
Research  and  NASA.  We  would  also  like  to  acknowledge 
useful  discussions  with  Drs.  H.  R.  Griem  and  M.  W.  C. 
Dharma-wardana . 


27 


V 


REFERENCES 


1.  H.  R.  Griera,  Proc.  7th  Int.  Conf.  on  Phenomena  in  Ionized  Gases, 

Belgrade,  Yugoslavia  (1965). 

2.  M.  Neiger  and  H.  R.  Griem,  Phys.  Rev.  A  14,  291  (1976). 

3.  J.  R.  Van  Zandt,  J.  C.  Adcock,  and  H.  R.  Griem,  Phys.  Rev.  A  lj», 

2126  (1976). 

4.  0.  Theimer  and  P.  Kepple,  Phys.  Rev.  A  1^,  957  (1970). 

5.  B.  F.  Rozsnyai,  J.  Quant.  Spectrosc.  Radiat.  Transfer  15^,  695 

(1975) }  Phys.  Rev.  A5,  1137  (1972). 

6.  S.  Skupsky,  Phys.  Rev.  A  2jL,  1316  (1980). 

7.  M.W.C.  Dharma-wardana,  private  communication  (1980). 

8.  M.W.C.  Dharma-wardana  and  R.  Taylor,  J.  Phys.  C:  Solid  St.  Phys. 

13,  XXX  (1980). 

9.  J.  R.  Greig,  H.  R.  Griem,  L.  A.  Jones,  and  T.  Oda,  Phys.  Rev.  Lett. 

24,  3  (1970). 

10.  D.  D.  Burgess  and  N.  J.  Peacock,  J.  Phys.  B4,  L94  (1971). 

11.  S.  Volonte ,  J.  Phys.  B8,  1170  (1975). 

12.  T.  L.  Pittman,  P.  Voigt,  and  D.  E.  Kellehex,  Phys.  Rev.  Lett.  45, 

723  (1980). 


28 


DISTRIBUTION  LIST 


David  D.  Burgess 
The  Blackett  Laboratory 
Imperial  College 
London  SW7  2BZ 
England 

Chander  P.  Bhalla 
Physics  Department 
Kansas  State  University 
Manhattan,  Kansas  66506 

J.  Cooper 
JILA 

University  of  Colorado 
Boulder,  Colorado  80302 

T.  R.  Carson 
University  Observatory 
St .  Andrews 
Fife,  Scotland 
UK  KY16  9LZ 

George  F.  Chapline 

Lawrence  Livermore  National  Laboratory 
P.0.  Box  808  L-71 

Livermore,  California  94550 

Raju  U.  Datla 

National  Bureau  of  Standards 
Building  221 

Washington,  D.C.  20234 
Richard  J.  Fortner 

Lawrence  Livermore  National  Laboratory 
P.0.  Box  808  L-401 

Livermore,  California  94550 

Aranon  Fisher 
Department  of  Physics 
University  of  California 
Irvine,  California  92717 

A.  H,  Gabriel 

Space  and  Astrophysics  Division 
Appleton  Laboratory 
Chilton,  Didcot,  Oxon. 

0X11  00X,  United  Kingdom 


Harold  C.  Graboske,  Jr. 

Lawrence  Livermore  National  Lab. 
P.  0.  Box  808  L-355 

Livermore,  California  94550 

Alan  Hauer 
M.S.  554 

Los  Alamos  National  Laboratory 
P.  0.  Box  1663 

Los  Alamos,  New  Mexico  87545 

Walter  F.  Heubner 
M.  S.  212 

Los  Alamos  National  Laboratory 
P.  0.  Box  1663 

Los  Alamos,  New  Mexico  87545 
Alan  L.  Hoffman 

Mathematical  Sciences  Northwest 

P.  0.  Box  1887 

Belevue,  Washington  98009 

Frank  E.  Irons 
Western  Mining  Corporation 
55  MacDonald  Street 
Kalgoorlie,  6430 
Western  Australia 

Lawrence  A.  Jones 
M.S.  455 

Los  Alamos  National  Laboratory 
P.  0.  Box  1663 

Los  Alamos,  New  Mexico  87545 

Elisabeth  Kallne 
Dept,  of  Physics 
University  of  Virginia 
Charlottesville,  Va.  22901 

Joseph  D.  Kilkenny 
The  Blackett  Laboratory 
Imperial  College 
London  SW7  2BZ 
England 

Kenneth  W.  Hill 
Plasma  Physics  Laboratory 
Princeton  University 
Princeton,  New  Jersey  08540 


Hans-Joach  I.  Kunze 
Inst.  Experimentalphysik 
Ruhr  Universitat 
Postfach  2148 
4630  Bochum 
West  Germany 

Richard  W.  Lee 
The  Blackett  Laboratory 
Imperial  College 
London  SW7  2BZ 
England 

James  G.  Lunney 
Department  of  Physics 
Queen's  University  of  Belfast 
Belfast,  Northern  Ireland 
United  Kingdom 

Eugene  J.  McGuire 
Org.  4211 

Sandla  Laboratories 
Albuquerque,  New  Mexico  87185 

Richard  M.  More 

Lawrence  Llvenuore  National  Laboratory 
P.  0.  Box  808  L-355 

Livermore,  California  94550 

Carl  Moser 

Centre  Europeande  Calcul  Atomique 
et  Moleculaire 
Batiment  506 
Universite  de  Paris  XI 
91405  Orsay  Cedex 
France 

Dennis  L.  Matthews 

Lawrence  Livermore  National  Laboratory 
P.  0.  Box  808  L-313 

Livermore,  California  94550 

Joseph  D.  Perez 

Lockheed  Palo  Alto  Research  Lab. 

3251  Hanover  Street  (5211) 

Palo  Alto,  California  94304 

G.  J.  Pert 

Dept,  of  Applied  Physics 
University  of  Hull 
Hull  HU6  7RX 
England,  United  Kingdom 


N.  J.  Peacock 
UKAEA 

Culham  University 
Abingdon,  Oxon 
OX  143DB,  United  Kingdom 

Jay  Pearlman 

Maxwell  Laboratories 

8835  Balboa  Avenue 

San  Diego,  California  92123 

A.  K.  Rajagopal 

Dept,  of  Physics  and  Astronomy 
Louisiana  State  University 
Baton  Rouge,  Louisiana  70803 

Balazs  F.  Rozsnyai 

Lawrence  Livermore  National  Laboratory 
P.  0.  Box  808  L-71 

Livermore,  California  94550 

Michael  Stroscic 
AF0SR/NP 

Bolling  Air  Force  Base 
Washington,  D.C.  203:2 

Stanley  Skupsky 

Laboratory  for  Laser  Energetics 
University  of  Rochester 
250  E.  River  Road 
Rochester,  New  York  14623 

David  Salzmann 

Soreq  Nuclear  Research  Center 

Yavne  70600 

Israel 

C.  Bruce  Tarter 

Lawrence  Livermore  National  Laboratory 
P.  0.  Box  808  L-71 

Livermore,  California  G4550 

Jon  C.  Weisheit 
Plasma  Physics  Laboratory 
Princeton  University 
P.  0.  Box  451 

Princeton,  New  Jersey  08544 
Wolfgang  Wiese 

National  Bureau  of  Standards 
Km.  A267 ,  Bldg.  221 
Washington,  D.C.  20234 


30 


David  M.  Woodall 
Dept,  of  Nuclear  Engineering 
University  of  New  Mexico 
Albuquerque,  New  Mexico  87131 

Barukh  Yaakobi 

Laboratory  for  Laser  Energetics 
University  of  Rochester 
250  E.  River  Road 
Rochester,  New  York  14623 

Zev  Zinamon 

the  Weismann  Institute  of  Science 
Rehovot ,  Israel 


Norman  Ro stoker 
Dept,  of  Physics 
University  of  California 
Irvine,  California  92717 


C .  F .  Hooper 
Dept,  of  Physics 
Univ.  of  Florida 
Gainesville,  Florida  32611 

John  D.  Hey 
Dept,  of  Physics 
Univ.  of  Cape  Town 
Rondebosch  7700 
South  Africa 

G.  Peach 

Dept,  of  Physics  &  Astronomy 

University  College 

Gower  Street 

London  WCIE  6BT 

England 

Naval  Research  Laboratory 
Code  4700  —  26  Copies 

Code  4720  —  50  Copies 


Alan  Kolb 

Maxwell  Laboratories 
8835  Balboa  Avenue 
San  Diego,  Calif.  92123 

Alan  Merts 

Los  Alamos  National  Scientific  Lab. 
P.  0.  Box  1663 

Los  Alamos,  New  Mexico  87545 


D.  E.  Kelleker 

National  Bureau  of  Standard 

Washington,  D.C.  20234 

H.  W.  Dravin 
EURAT0M-CEA 

Dept,  of  Plasma  Physics 
Center  for  Nuclear  Study 
F-92260  Fonteney-aux-Roses 
France 


31 


